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The structure and electronic order at the cleaved (001) surfaces of the newly- discovered pnictide 
superconductors BaFe2-xCoxAs2 with x ranging from to 0.32 are systematically investigated by 
scanning tunneling microscopy. A a/2 x a/2 surface structure is revealed for all the compounds, and 
is identified to be Ba layer with half Ba atoms lifted-off by combination with theoretical simulation. 
A universal short-range charge order is observed at this a/2 x a/2 surface associated with an energy 
gap of about 30 meV for all the compounds. 

PACS numbers: 73.20.-r, 68,37.Ef, 74.25.Jb 



The newly discovered high-Tc superconductors based 
on iron pnictides have ignited intense research interests, 
following the cuprate boom [l|, 0, i, i, i, 0, 0- Al- 
though progresses have been made to understand their 
structural, electronic, and magnetic properties, the na- 
ture of superconductivity is yet a mystery [8]. Among 
the experimental techniques used to investigate the 
high-Tc superconductors, surface-sensitive probes such 
as angle-resolved photoemission spectroscopy (ARPES) 
and scanning tunneling microscope (STM) have played 
a critical role, due to their unrivaled energy and spa- 
tial/momentum resolution. However the surface might 
differ from the bulk structurally and electronically, due 
to the symmetry-breaking and/or reconstruction. So 
thorough understanding of the surface structure and 
electronic properties at atomic scale is crucial to inter- 
pret correctly the experimental data based on surface- 
sensitive detections. On the other hand, the discrepant 
correlation between lattice, charge, and spin at the sur- 
face may lead to novel two-dimensional collective order, 
which would open a new door to the nanoelectronics ap- 
plications based on correlated materials. 

STM has been used to investigate the 122-type pnic- 
tide compounds, for which high-quality single crystals 
can be srown 0,0. Stripe |9|,[M[ni or square-lattice 
[12I , [13I , [l3 | structures were observed at the surfaces 
of some certain compounds. However, the surface ter- 
mination of the 122-type pnictides is still in debate, 
and the evolution of the geometric and electronic struc- 
tures as a function of doping level has never been ex- 
plored at atomic scale. In this letter, we report sys- 
tematic STM study of electron-doped 122 compounds 
(BaFe2-a:Coa^As2) with nominal doping composition of x 
= 0, 0.08, 0.17, 0.20, 0.25, 0.32. A universal ^2x^2 sur- 
face structure is observed for all the doping levels. With 
the aid of theoretical simulation, this surface structure 
is determined to be well-ordered Ba plane with half Ba 
atoms lifted-off. Strikingly, a short-range charge order is 



characterized on such x ^/2 surface for all the com- 
pounds accompanied by an energy gap ranging from 24 
to 40 meV. 

The growth of pnictide single crystals and their mag- 
netic/transport characterization have been described 
elsewhere [15|. The crystals were cleaved in-situ at ~120 
K in ultra-high vacuum environment with pressure lower 
than 2xl0~^^ mbar, and immediately transferred into 
the STM stage, which were already cooled to 5 K. 

Two typical structures have been observed at the 
cleaved surfaces for all the investigated samples. The 
first is the widely observed stripe structure (not shown) 
[1, The other type of surface is atomically flat, 

as shown in Fig. 1(a). The typical images with atomic 
resolution, for instance, for x = and x = 0.20, are 
shown in Fig. 1(b) and (c), which have square lattice 
with lattice constant of ~5.6 A . This observed universal 
square-like structure is quite striking, since the structure 
is orthorhombic for x < 0.20, while tetragonal for x > 
0.20 at 5 K The typical STM images with atomic 

resolution for x = (parent) and x = 0.20 at 5 K are 
shown in Fig. 1(b) and (c), and both show 5.6 A x 5.6 
A lattice. Scanning tunneling spectroscopy (STS) mea- 
surement was also performed on the cleaved surfaces with 
square-like lattice for all the investigated samples at 5 K, 
and the typical dl/dV spectroscopies taken away from 
defects are shown in Fig. 1(g). A protrusion feature at 
about -0.2 V is observed in the STSs for all the samples, 
which could originate from collective ordering or band- 
structure effect. 

In the BaFe2-icCoa:As2 compounds consisting of FeAs 
layers separated by a single Ba layer, the bond between 
Fe and As atoms is quite strong, so the cleavage is likely 
to happen between the FeAs layer and Ba layer, or in the 
Ba layer, leaving half of the Ba atoms on each exposed 
surfaces [9]. The latter seems more preferable energet- 
ically since this cleavage can still maintain charge neu- 
trality in the system. In order to determine the universal 
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surface structure, we did theoretical calculation with half 
layer of Ba atoms at the surface, which form \/2 x a/2 
structure (in the tetragonal notation). Our calculations 
are carried out by employing the Vienna ab initio simula- 
tion package (VASP) using the general gradient approxi- 
mation (GGA) with PBE functional and plane wave basis 
sets [13, mi, [lH, Eoj . A supercell geometry with 9 atomic 
layers (4 FeAs layers and 5 Ba layers) and a vacuum layer 
about 30 A thick are used to simulate the surface. The 
plane wave cutoff is set to be 600 eV. For geometry opti- 
mization, we use a 4x4x1 Monkhorst grid [21] to sample 
the Brillouin zone, while a denser 6x6x1 grid is used 
for density of states calculations. During the structural 
optimizations, we fix the bottom 3 (2 Ba and 1 FeAs) lay- 
ers in the bulk configuration and allow all other atoms 
in the supercell to move until all forces vanished within 
0.02 eV/A . 

The simulated STM images for the parent with or- 
thorhombic and tetragonal structures are displayed in 
Fig. 1(d), and both show square-like pattern. No es- 
sential difference can be resolved from our simulation, 
although the symmetries of the two structures diverge. It 
is reasonable considering that the distortion of the low- 
T orthormhombic structure from the high-T tetragonal 
structure in BaFe2As2 is quite minor [l,[l^. The simu- 
lated STM image for x = 0.25 (tetragonal structure) is 
shown in Fig. 1(e), which is similar to that of the parent 
compound. Orthorhombic distortion also has little ef- 
fect on the simulated images. So the simulated STM 
images with half Ba layer termination well reproduce 
the experimental data for both undoped and doped com- 
pounds. STM simulation for the parent compound with 
FeAs layer as the surface termination was also performed, 
and only shows a square-like lattice of 4 A x 4 A (or 1 x 
1 in the tetragonal notation). This is inconsistent with 
the observed STM images. 

We also calculated the local density of states (LDOS) 
of the parent for different layers, and there is indeed a 
feature arising around -0.2 eV for the first Ba and FeAs 
layer, which is consistent with the STS results. The 
agreement between experimental STM/STS results and 
theoretic simulations strongly suggests that the cleaved 
surface for the series of BaFe2-a:CoajAs2 compounds with 
X ranging from to 0.32 is terminated by half layer of Ba 
atoms. The half Ba atoms form well-ordered a/2 x a/2 
structure, which is insensitive to the orthorhombic struc- 
tural distortion. 

STS in narrower range around Fermi level was mea- 
sured on the a/2 x a/2 surface for all the investigated 
samples, and the typical spatially resolved dl/dV spec- 
troscopies obtained away from defects for x = 0, 0.20, and 
0.32 are shown in Fig. 2. The dl/dV curves on all the 
samples show similar shapes, and surprisingly a weak gap 
structure is observed for all the samples. The gap size 2A 
is similar for all the samples, ranging from 24 to 40 meV. 
This gap size is much larger than those previously ob- 
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FIG. 1: (a) Large-scale STM image of the BaFe1.8Coo.2As2 
surface with ^2 x a/2 structure (Vs = 0.5 V, I = 0.5 nA). 
(b) and (c) are STM images showing the \/2 x >/2 lattice 
with atomic resolution on the surfaces of BaFe2As2 (Vs = 
30 meV, I = 3 nA) and BaFe1.8Coo.2As2 (V^ = 2 meV, I = 
20 nA), respectively, (d) and (e) are the simulated empty- 
state STM images with half Ba atoms on top of BaFe2As2 
and BaFei.75Coo.25As2, respectively, (f) The side view and 
top view of the calculated structure of BaFei.75Coo.25As2. 
The calculated structure is orthorhombic in (d), and tetrag- 
onal in (e) and the inset of (d). (g) dl/dV spectroscopies for 
BaFe2-a;CoxAs2 with x = 0, 0.08, 0.17, 0.20, and 0.25, respec- 
tively, (h) Calculated LDOS for different layers in BaFe2As2. 

served by STM (-12 meV) EEI and ARPES (13 and 
10 meV for different Fermi surfaces) [22[ on the striped 
surface, which were attributed to the isotropic s-wave su- 
perconducting gap. The universality of the observed gap 
feature on the x ^/2 surfaces strongly indicates that 
they share the same origin. Superconducting gap due to 
bulk superconductivity is excluded since the gap was ob- 
served in the parent compound. Neither can the band 
structure of the \/2 x \/2 surface termination give rise 
to such a gap from our calculation. It also can not be 
attributed to the bulk SDW order, as SDW disappears 
when X > 0.20 [15j]. It could be the manifestation of 
some surface collective order induced by the new 
surface termination, such as surface superconductivity, 
surface spin order or surface charge order. Occasionally 
gap feature becomes very weak at some random locations 
(blue curves in Fig. 2). The observation of universal gap 
feature indicates that the electronic correlation on the 
surfaces depart significantly from the bulk, and that the 
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FIG. 2: Dl/dV spectroscopies measured along line trajecto- 
ries on the BaFe2-xCoxAs2 compounds with x = (a), 0.20 
(b), and 0.32 (c). The lengths of the trajectories are 31 nm, 
32 nm, and 28 nm, respectively. The purple and blue curves 
represent spectroscopies with and without gap feature, re- 
spectively. 
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prediction of bulk properties based on surface-sensitive 
probes should be very careful. 

In order to explore possible charge ordering on the 
surface, dl/dV mapping simultaneously with topography 
imaging was performed on all the investigated samples 
at different bias voltage. Typical STM images and the 
corresponding dl/dV maps for x = 0.17 are shown in 
Fig. 3(a) and (f) with sample bias of 0.02 V and -0.02 
V, respectively. The >/2 x >/2 lattice with atomic reso- 
lution is clearly resolved in both topography and dl/dV 
images at such low tunneling-j unction resistance. Fur- 
thermore, inhomogeneous corrugations at nanoscale are 
observed, superimposed on the square lattice. This is 
also evidenced in Fig. 1(b) and (c). If the tunneling cur- 
rent is lowered while the sample bias is unchanged, the 
atomic resolution becomes weaker as shown in Fig. 3(b) 
and (g), due to the increased distance between the tip 
and the substrate, while the nanoscale corrugations are 
still preserved. When the bias voltage is swept, the spa- 
tial inhomogeneity pattern in the same area also varies 
correspondingly, and even reverses at some locations as 
clearly shown in Fig. 3. So the nanoscale inhomogene- 
ity is unlikely to originate from topographic effect, and 
should be due to electronic effect. 

The larger-scale STM images and corresponding dl/dV 
maps for different doping concentrations are shown in 
Fig. 4. It is clear that the inhomogeneous charge or- 
der is a universal phenomenon, which can be observed 
for all the compounds, including the parent compound, 
although it is relatively weak. The inhomogeneity pat- 
tern does not have spatial long-range order as indicated 
by the fast Fourier transformation (not shown), nor ob- 
vious preferable orientation. This short-range charge 
order resembles the spatial charge fluctuation in YSi2 
nanowires [l^ . The evolution of the charge inhomogene- 
ity as a function of the sample bias is also investigated. 
It is interesting that the spatial pattern and character- 



FIG. 3: A series of topography images (upper) and corre- 
sponding dl/dV maps (lower) measured at the same area with 
the size of 6.25 nm x 6.25 nm when varying the tunneling pa- 
rameters: (a) Vs = 0.02 V, I = 2 nA; (b) = 0.02 V, I = 
0.5 nA; (c) V. = 0.05 V, I = 0.5 nA; (d) V. = 0.1 V, I = 0.5 
nA; (e) = 0.4 V, I = 0.5 nA; (f) V. = -0.02 V, I = 1.5 nA; 
(g) Vs = -0.02 V, I = 0.5 nA; (h) = -0.05 V, I = 0.5 nA; 
(i) Vs = -0.1 V, I = 0.5 nA; (j) = -0.4 V, I = 0.5 nA. 



istic length of the charge inhomogeneity vary when the 
voltage is changed. We did autocorrelation calculation 
for the dl/dV images for quantitative analysis. The cor- 
relation lengths can be extracted from the radical line 
profiles, which characterize the mean distance between 
the neighboring maxima/mimima [ij]. The correlation 
length depending on the doping level (x = 0.17, 0.20, 
and 0.25) and energy is shown in Fig. 4(i). Surprisingly, 
the correlation length at a given energy and its evolu- 
tion with energy is quite similar for different doping lev- 
els, although they show different order and/or ordering 
strength in the bulk. The averaged correlation length is 
about 3.0 nm. Since the charge inhomogeneity and en- 
ergy gap are universally observed on the \pl x \pl surface 
for all the investigated compounds with a wide doping 
range, it is natural to think that they are connected, i.e., 
the energy gap reflects the characteristic energy of the 
short-range charge order. 

For cuprates, spatially periodic charge modulations 
have been seen within the energy gap due to quasiparticle 
interference oscillations in the superconducting state 1251 
or electronic modulations in the pseudogap state 2^, 27| . 
The periodicity is energy dependent for the former and 
independent for the latter. We acquired topography and 
dl/dV images with sample bias of -0.01 V, within the en- 
ergy gap as displayed in Fig. 4(j). The observed charge 
pattern does not show any regular features. This might 
reflect the different pairing symmetry of pnictides (ex- 
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FIG. 4: Large-scale topography (upper) and dl/dV (lower) 
images of BaFe2-xCoxAs2 for x = (a, b), 0.17 (c, d), 0.20 
(e, f), 0.25 (g, h). The sample bias is 0.1 V for (a), (c), (e), 
(g) and -0.1 V for (b), (d), (f), (h). The tunneling current is 
0.5 nA for all the images. The images for the same doping 
concentration are in the same area, (i) Correlation length 
derived from the dl/dV images as a function of sample bias 
for X = 0.17, 0.20 and 0.25. (j) Topography and dl/dV images 
for X = 0.17 at V, of -0.01 V. 



tended s-wave probably [28]) from cuprates (d-wave). 

DFT calculation failed to reproduce such charge inho- 
mogeneity in both the parent and doped compounds as 
shown in Fig. 1(b) and (c). It has been predicted the- 
oretically that electron interactions in iron pnictides are 
sufficiently strong to produce incipient Mott physics [2^ . 
The charge order might come from the enhanced elec- 
tron correlation effect at the symmetry-broken \/2 x \/2 
surface, though more efforts are needed to identify the 
microscopic origin. 

In conclusion, a >/2 x \/2 surface structure terminated 
by half layer of Ba atoms is identified at the surfaces 



of BaFe2-a:Coa^As2 compound with x varying from to 
0.32, which is insensitive to the orthorhombic distortion. 
A universal short-range charge order accompanied by an 
energy gap of around 30 meV appears at this x ^/2 sur- 
face for all the investigated Co doping levels. This univer- 
sal charge order suggests relatively strong electron cor- 
relation in the iron pnictides, at least on the symmetry- 
broken surfaces. It also sheds light on the exploration 
of novel cooperative effects of the correlated materials in 
reduced dimensions. 
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